function phi=phi_n(x,n)
x_max = max(x);
phase = rand()*pi;
% phase = 0.0;
phi = sqrt(2/x_max)*sin(2*n*pi*x/x_max+phase);
